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ABSTRACT 

We present BIMA observations of a 2' field in the northeastern spiral arm of M31. In this region 
we find six giant molecular clouds that have a mean diameter of 57±13 pc, a mean velocity width of 
6.5±1.2 km s , and a mean molecular mass of 3.0 ± 1.6 x lO 5 A/ . The peak brightness tempera- 
ture of these clouds rang es from 1.6-4.2 K. We compare these clouds to clouds in M33 observed by 
iWilson fc ScmoTiel (fl99(f using the OVRO millimeter a rray, and some cloud complexes in the Milky 
Way observed bv iDame. Hartmann fc Thaddeusl (|2001l ) using the CfA 1.2m telescope. In order to 
properly compare the single dish data to the spatially filtered interferometric data, we project several 
well-known Milky Way complexes to the distance of Andromeda and simulate their observation with 
the BIMA interferometer. We compare the simulated Milky Way clouds with the M31 and M33 data 
using the same cloud identification and analysis technique and find no significant differences in the 
cloud properties in all three galaxies. Thus we conclude that previous claims of differences in the 
molecular cloud properties between these galaxies may have been due to differences in the choice of 
cloud identification techniques. With the upcoming CARMA array, individual molecular clouds may 
be studied in a variety of nearby galaxies. With ALMA, comprehensive GMC studies will be feasible 
at least as far as the Virgo cluster. With these data, comparative studies of molecular clouds across 
galactic disks of all types and between different galaxy disks will be possible. Our results emphasize 
that interferometric observations combined with the use of a consistent cloud identification and anal- 
ysis technique will be essential for such forthcoming studies that will compare GMCs in the Local 
Group galaxies to galaxies in the Virgo cluster. 

Subject headings: galaxies: Local Group (M31, M33, Milky Way) — galaxies: ISM — radio lines: ISM 



1. INTRODUCTION 

In the inner disk of the Milky Way, the dominant com- 
ponent of the interstellar medium is molec ular gas which 
is dis t ributed i n disc rete cloud complexes (|Scoville et al.l 
Il987t IScovilld I1990D . Nearly all star for mation in our 
Galaxy is assoc iated with these clouds (|Scoville et al.l 
ll987tlBlitziri993h . Therefore, understanding the distribu- 
tion and properties of molecular clouds is a prerequisite 
for understanding galaxy evolution. Studies of the molec- 
ular gas emission in the Milky Way show that while the 
cloud mass spectrum extends over several orders of mag- 
nitude, a majority of the molecular gas mass is contained 
in large, massive cloud complexes, which are typically 
^40 pc in diameter and have masses on the order of 1 
x 10 5 M Q ; these complexes are usually referred to as gi- 
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ant molecular clouds (GMCs) 1 and are considered to be 
the basic organizationa l unit of the molecular interstellar 
medi um (lScovilldll990t lCombedll99l1: lYoung fc Scc^fill 

EMU). 

Compared to the Milky Way, the properties of GMCs 
in external galaxies are not well known. The primary rea- 
son for this is a lack of high resolution, high sensitivity 
data which are necessary for spatially resolving individ- 
ual complexes. Only in the Magellanic Clou ds can single 
dish t elescopes resolve extragalactic GMC s (jlsrael et al.l 
119931: iRubio et al.lll993t iFukui et al.ll2001h . Even in the 
nearest spirals (e.g., M31 or M33), the angular size of 
a typical 40 pc GMC is ~12", approximately one half 
the resolution of the IRAM 30m single dish telescope at 
2.6mm, the wavelength of the CO (J=l-0) line emission. 
For these galaxies, interferometric observations are nec- 
essary to resolve GMCs. 

GMCs have been detected using millimeter in- 
terferometers in the two nearest spiral galaxies, 



1 In this paper, we will use the words giant molecular clouds, 
GMCs and cloud complexes interchangeably. 
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19871 IWilson fc Rudolphl 119931 lLoinard fc Ailed 1 19911 
Rosolowskvl 12007ft . In M33. IWilson fe Scovilld (<1990T ) 



observed 19 different fields and identified over 30 GMCs 
but only 9 GMCs were suitable for a study of cloud 
properties. The number of GMCs detected in most spiral 
galaxies has been less than a dozen, presumably due to 
the small coverage area. This is beginning to change. 
Recently a survey of M33 was undertaken by the BIMA 
array. The angular resolution was c oarse (~13") but 
adequ ate to iden tify nearly 150 GMCs (lEngargiola et al.1 
120031 see also IRosolowskv et al.l 120071 ). With the 
upcoming CARMA array, ground-breaking studies like 
the M33 survey are likely to become routine for external 
galaxies. 

It is important to note that previous studies of GMCs 
in external galaxies have used different instruments with 
varying angular and linear resolutions and varying sensi- 
tivity. Moreover, these studies have used different cloud- 
identification methods to define GMCs and determine 
their properties. Yet, in almost all of these studies, the 
observations have detected discrete molecular features 
similar to Galactic GMCs. Particularly notable is the 
consistent linewidth-size relationship in most of these ob- 
servations. 

However, several studies have noted differences be- 
tween Galactic GMCs and those found in some irreg- 
ular and dwarf ga laxies. In the i rregular galaxies NGC 
6822 (jWilsonlll994ft and the SMC (jRubio et al.lll993h . the 
clouds are slightly less massive, and perhaps smaller than 
clouds seen in the Milky Way. In the LMC. ICohen et aTl 
(1988) find that for a given size and line width, clouds 
are about six times fainter in CO than comparable clouds 
in the Milky Way. The same is true of the dwarf NGC 
1569 where the CO-to-H2 conversion factor is reported to 
be six times larger than the Galactic value (jTavlor et al.l 
1999). In IC10, the cloud sizes and masses are similar to 
the Milky Way clouds, but the CO-to-H 2 conversion fac- 
tor may be twice as high as that seen in t he Galaxy if the 
cloud s are assumed to be self-gravitating (jWilson fc Reidl 
|1991|) . These studies suggest that there are noticeable 
differences between the Galactic GMCs and those found 
in irregular or dwarf galaxies, that may be attributed 
to the different host environments (e.g., metallicity) of 
these galaxies. 

GMCs in the three nearest spiral galaxies, M31, M33 
and the Milky Way can be observed at sufficiently high 
resolution (~20 pc) to remove ambiguities from coarse 
resolution. Are the GMC prop erties in these spira l galax- 
ies fundamentally different? IWilson fc Scovillg (|1990f ) 
claim that no GMC larger than 4x 10 5 M is seen in M33, 
and unlike the Milky Way, 50% of the molecular mass in 
M33 resides in small clouds with masses less than 8xl0 4 
Mq, based on a comparison of interferometric and single 
dish flux from the CO(J=1-0) line. On the other hand, 
analysis of the 13 CO/ 12 CO line ratio for a number of M33 
clouds suggests that the total mass measured from the 
12 CO (J=l-0) line may be overestimated and the actual 
conversion factor from the CO to H2 mass may be lower 
in the diffuse gas; in that case, the amount of mass in 
GMC s would constitute most of the molecular mass in 
M33 ([Wilson fc Walkerl Il99i) With coarser resolution 
data, Engargi ola et al.1 (|2003ft find a steep GMC mass 



spectrum in M33 and report a characteristic mass of 7 x 
10 4 M Q . They conclude that molecular clouds in M33 are 
formed rapidly from the atomic gas and are short lived 
entities. 

On the other hand, in M 31 IVogel et al.l (|1987l ) 
and IWilson fc Rudolphl (|1993l ) note that the four 
GMCs they d etect are similar to Galactic GMCs. 
lLoinard et alj (|1999h smooth a section of the Ca- 
rina arm in the Milky Way to the resolution of sin- 
gle dish M31 maps and conclude that there are gen- 
eral similarities between the distribution of molecu- 
lar gas in M31 and the Milky Way. Similar results 
were also found by iHever. Dame, fc Thaddeusl (|2000f ) 
who compared the molecular gas in M31 and the 
Milky Way using single dish o bserva tions; in their 
studv lHever. Dame, fc Thaddeusl (|2000ft find a number 
of similarities in the molecular gas proeprties between 
M31 and the Milky Way such as the amount of molecu- 
lar gas at radii larger than 8 kpc, large arm to interarm 
contrasts, and similar surface bri ghtness, line widths and 
spacings of GMCs in spiral arms. iRosolowskvl (|2007f) has 
recently completed a large survey of clouds in M31. He 
also finds that the GMCs in M31 are similar to those 
found in the Milky Way. These results depend on com- 
parison of properties derived from datasets which are ob- 
tained with a different instrument and analyzed with dif- 
ferent methods of cloud identification. The Milky Way 
data, for example, come from single dish observations 
with a signal to noise that is at least 2-4 x higher than 
in the M33 and M31 datasets. The datasets also dif- 
fer in that single-dish telescopes are inherently different 
than the aperture synthesis maps: the interferometer will 
miss smooth, extended emission whereas the single-dish 
data may suffer from beam dilution. The Milky Way 
survey also suffers from varying linear resolution. Fi- 
nally, the methods for identifying clouds in the different 
Milky Way studies vary. T hese studies (e.g.. lDame et al.l 
1986; Sanders et al.lll985ft have usually projected a typi- 
cal (l,b,v) cube on one of its axes and used a integrated 
intensity contour to identify discrete features they call 
clouds or complexes. The sizes were measured from the 
total are a of the cloud, as suming that the cloud was 
spherical (jDame et al.l ll986). or by measuring the chord 
across the velocity ce ntroid in a position- velocity diagram 
dSanders et al.lfl985h . In M33, (IWilson fc Scovilld fl990l) 
used a two-tiered selection criterion which required that 
a cloud be 3er or more in 2 adjacent channels and that 
its summed flux be at least 3a above the noise. 

The lack of consistency in these methods and differ- 
ences in observing techniques are worrisome. A proper 
comparative study of GMCs in these galaxies requires 
that the differences described above be eliminated. This 
is particularly important in view of future telescopes such 
as CARMA and ALMA which will have the ability to 
mosaic large regions of galaxies like M31 and M33. The 
large scale surveys will have the detail necessary for in- 
depth comparisons of the molecular ISM in Local Group 
and more distant galaxies. These future studies will ad- 
dress the fundamental question of how the molecular ISM 
varies from galaxy to galaxy. Knowing this is critical to 
understanding star formation and galaxy evolution. 

In this paper, we make an attempt to do a fair compar- 
ison of the molecular emission in the three nearest spiral 
galaxies, i.e., M31, Milky Way and M33. We compare 
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the M31 complexes to the simulated Milky Way com- 
plexes and previous surveys, using a cloud- identification 
technique in which we identify complexes using a single 
integrated intensity contour (<j3]), and discuss how pre- 
vious conclusions about M33 would change if the same 
technique of cloud identification were used ( §4.1jl . An im- 
portant addition to this analysis would be the addition 
of single dish data for the M31 and M33 clouds which 
would allow for an even more comprehensive compari- 
son. Of course, one would still need to project the Milky 
Way clouds to the distance of Andromeda to compare the 
clouds at the same spatial resolution. Our method is in- 
tended to lay the groundwork for future studies that will 
compare molecular clouds across galactic environments 
and across different galaxy types with telescopes such as 
CARMA, SMA and ALMA. 

In ^4.21 we describe two automated c loud identification 
algori thms (Gaussclump developed bv lStutzki fc Glisten 
( 1990D . and Clumpfind developed by I Williams et al. 
(1994)), and their advantages and disadvantages. While 
we do not advocate using these algorithms to describe the 
characteristics of GMC populations in a galaxy, these al- 
gorithms may be used to compare the M31, M33 and 
the simulated Milky Way data. We describe the results 
of this experiment in £14 . 31 and discuss our overall conclu- 
sions in fj5] Preliminary results from this work have been 
published in lSheth et all <|2000D . 

2. OBSERVATIONS 

2.1. M31 BIMA Observations 

We observed CO(J=1-0) emission in the northeast- 
ern spiral arm of M31 with the BIMA (Berkeley-Illinois- 
Maryland Association) 2 array. Using a 7-pointing hexag- 
onal grid (one central pointing surrounded by a ring of 
six pointings, each separated by 53" from its two neigh- 
bors in the ring and the center), we mosaicked a region 
kn own to have strong si ngle-dish CO emission (point 4,0 
in (|Rvden fc St ark 1986)). The mosaicked fields were ob- 
served sequentially with an integration time of 1 minute 
per field for 35 minutes, at the end of which, the calibra- 
tor was observed for 3 minutes. The total data presented 
here was collected in five separate tracks for a total in- 
tegration time of 18.2 hours, or 2.6 hrs per field . The 
fields were combined using MI RIAD (ISau t et al.Hl995h 
with an algorithm developed bv lSault et al.1 ( 19961 ). Note 
that the final maps generated from this algorithm have 
variable gain and sensitivity. The noise increases mono- 
tonically from the map center towards the edge. The 
gain is constant over a fairly large region (in this case, 
the central ^180") and then drops sharply to the edge. 
For our analysis, we only considered the region where the 
gain was higher than 0.4. Further details of the observa- 
tions are given in Table [TJ 

2.2. M33 and Milky Way Data 

The M33 survey was conducted bv IWilson fc Scovilld 
|l990h with the OVRO (Owens Valley Radio Observa- 
tory) milli meter array. The data r eduction process is 
described in lWilson fc Scovilld (j 19901) . 

2 The BIMA Array was partially funded by grant AST-9981289 
from the National Science Foundation 

3 The OVRO Array was partially funded by grant AST-9981546 
from the National Science Foundation 
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Fig. 1. — Velocity-integrated, primary beam corrected, intensity 
map of the M31 field shown in logarithmic contours of 1.5 n x 2 K 
km/s, where n=-l, 1,2,3,... The grey scale image is a DSS J-band 
image in reverse color, so that white color represents obscuring 
dust. The six complexes used for comparison to the Milky Way 
are identified with the boxes labelled A-F. The outer circle is the 
gain=0.4 contour. 

The Milky Way observ ations and reductions are from 
the 1.2m CfA survey (|Dame. Hartmann fc Thaddeusl 
2001). The different complexes studied in the present 
work have been previously pub lished as fo llows: Gem 
OBI dStacv fc Thaddeusl fl99lT). W3CO (iDieel et al. 
1996), Cas A (|U neerechts. Umb anhowa r. fc Thaddeusl 
| 2000li. and the Carina arm d Grabelskv et all 119871: 
iBronfman et alj Il989t ) . Table [H lists the datasets and 
instruments used in this study. 

3. RESULTS: M31 CLOUDS 

In Figure [TJ we show the velocity- integrated CO(J=l- 
0) intensity map in contours overlaid on a grey-scale opti- 
cal 4 image. Our field of view encompasses two dust lanes 
which can be seen trending diagonally across the field 
from NE to SW. The majority of the CO emission lies in 
these two dust lanes, which are separated by a distance 
of 2 kpc (assuming an inclination, i=77°). There are six 
distinct molecular cloud complexes, which are labeled A- 
F. The complex labeled E lies between the two dust lanes 
along a dust spur which runs between the two dust lanes. 
We rotated the data cube, collapsed one spatial dimen- 
sion and generated RA-velocity and DEC-velocity maps 
for each complex. We checked these maps to verify that 
each complex is a well-defined and separate entity. None 
of the complexes have emission separated by more than 
a few km s _1 except complex C. 

In Figure^ we show channel maps in which CO (J=l- 
0) emission is detected. These maps highlight the com- 
plex structure present in these clouds. Complex A, to- 
wards the northern edge of our field of view, is in the dust 
lane adjacent to a bright HII region. Further south in the 
dust lane is complex B which has considerable structure; 
it is made up of at least two or three different clumps 
which the automated algorithms separate into different 

4 Digiti zed Sky Survey image of M31 obtained from the Sky view 
database (McGlynn e t al.H1996h 
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Fig. 2. — Channel maps showing the CO (J=l-0) emission in the 
northeastern spiral arm of M31. These maps are not primary beam 
corrected. The contours are logarithmic, 1.5 n x 0.22 Jy beam -1 . 
The boxes show the complexes A-F. The channel maps are 2.6 km 
s _1 wide starting at -142 km s -1 . The outermost circular contour 
is the gain=0.4 contour. 

clumps. Complex C is the weakest complex in our maps. 
It is patchy and extends over a large region, and over 10 
km s _1 in velocity. The spectrum for this complex shows 
that it may be a blend of several clouds. Therefore we do 
not calculate the size or velocity width for this complex. 
Complex D may be considered a prototypical GMC; it 
is well-isolated, i.e. unblended and has nice symmetry. 
This complex is also rather weak but we are fortunate 
that it is near the center of our mosaicked field where 
the noise is the lowest in these maps. Complex E lies 
along the dust spur connecting the two dust lanes and 
also has considerable structure. There are at least two or 
three peaks of emission in this complex. Complex F is an 
elongated complex towards the western edge of our field 
of view; this is the brightest complex in our field of view. 
We also show the boxes used to extract the spectrum and 
sum the emission for each complex. From these, we cal- 
culated the virial and molecular mass for each complex. 
The properties thus derived are listed in Table [5] The 
method used for identifying these clouds is discussed in 
the next section. 

4. COMPARING CLOUDS WITH DIFFERENT CLOUD 
IDENTIFICATION ALGORITHMS 

4.1. The Integrated Intensity Contour Method 

We used a specified integrated intensity contour (2 K 
km s _1 ) to define a complex. Is this justified, espe- 
cially when many of these complexes clearly have sub- 
structure? From our experience with previous Milky 
Way studies, we believe that any cloud complex can 
always be divided and sub-divided into smaller and 
smaller clouds because the molecular medium is ex- 
tremely clumpy with a low volume filling factor (see e.g., 
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Fig. 3. — The size-linewidth relationship is shown for the M31 
complexes (stars), simulated Milky Way GMCs (squares, circles, 
circles with a plus in the mi ddle) and previous ly analyzed Milky 
Way complexes (triangles) by Dame et al. ( 1986). The typical error 
in measurements of the velocity width and radius are shown with 
the cross at the mean value of the M31 clouds in the survey. All 
the data points fall on the linewidth-size relationship defined by 
the Dame et al. (1986) Milky Way data. The data, whether from 
a single dish telescope or an inteferometer, gives equivalent results 
when the observations sample the same spatial scales and have 
comparable sensitivities in different galaxies. 



iDame. Hartmann fc Thaddeusl [200lT ). The higher the 
resolution, the more likely it is that a given structure may 
be divided into smaller components. Having said that, 
however, a complex has definite boundaries. Therefore 
the integrated intensity contour in a collapsed cube (i.e., 
collapsed along one axis) seems to be a reasonable way 
of identifying a molecular cloud complex and evaluating 
its properties. Thi s has been applied to GMC studies 
in the Milky Way (|Sanders et al.lll985l . e.g.,). Thus we 
can compare M31 clouds directly to the previous Milky 
Way studies. As a check on the method, we applied 
the method to a few Milky Way complexes (Gem OBI, 
W3CO, Cas A) and to a few complexes in M33. We 
projected the Milky Way clouds to the distance of M31 
and simulated their interferometric observation using the 
same uv tracks actually obtained for the BIMA observa- 
tions of M31. The details of these simulations are in the 
Appendix. Using this approach each of the three Milky 
Way complexes was decomposed into two clouds. We 
refer to the Milky Way complexes projected to the dis- 
tance of M31 and observed with the interferometer as 
"simulated" clouds for the remainder of this paper; the 
reader should keep in mind that these clouds are not ar- 
tificial but are actual Milky Way clouds which are being 
"observed" as if they were in Andromeda and were ob- 
served with BIMA. The properties of the complexes thus 
defined are listed in Table 2J 

When we compare the properties of the complexes 
we observe in M3 1 to those observed previously in th e 
Milky Way (e.g.. IDame et al.lll986l ; ISanders et al.lll985h . 
we find that the mean properties of the complexes in 
the two galaxies a re not significantly different. Whereas 
IDame et al.l ()1986f ) found complexes which ranging from 
24-262 pc in diameter with velocity widths of 4-18 km 
s -1 , the complexes in M31 range from 40-75 pc and 4.7- 
7.8 km s _1 . In Figure [3 we place our M31 com plexes, 
the Milky Way clouds from IDame et all (|1986l ). and a 
set of sim ulated Milky Way c louds o n the size-linewidth 
plot from IDame et alj (|1986h . The IDame et al.l (|1986l ) 
GMCs are particularly useful for this figure because they 
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have a larger dynamic range in size and linewidth than 
the other sets of clouds. Within the errors, the M31 
complexes fall along the same relationship found in the 
much more extensive Milky Way survey. It is also worth 
noting that the Milky Way GMCs simulated at the dis- 
tance of Andromeda and observed with an interferometer 
show linewidths and sizes that are also consistent with 
the standard relationship. The main point of this figure 
is that data for different galaxies whether observed with 
a single dish telescope or an interferometer can be com- 
pared equivalently so long as the observations have com- 
parable spatial resolution (i.e., the observations sample 
the same spatial scales in different galaxies) and sensi- 
tivity, and the data are analyzed with the same cloud 
identification technique (s). 

For identifying clouds in M33, iWilson fc Scoviliel 
(fl99l did not use the same cloud identification method 
as described above. In their method, if a cloud could be 
resolved by eye in a channel and separated spatially or in 
velocity from an adjacent cloud, then the cloud was con- 
sidered to be a distinct cloud. The effect of this technique 
is similar to the automated algorithms like Clumpfind 5 
and Gaussclump which we discuss in more de tail later 
in SU Basically, the I Wilson fc Scoviliel (jl990l ) method 
breaks apart complexes whenever possible. Thus it was 
difficult to deduce from their analysis whether the M33 
clouds are indeed a different population than the Milky 
Way clouds. 

We re-analyzed Figure 2 from IWilson fc Scovi 11c (1990) 
using the integrated intensity contour method and found 
that several of the M33 clouds merged into larger com- 
plexes. The main reason for that is that the lowest c on- 
tour shown in Figure 2a of IWilson fc Scoviliel (|1990l ) is 
7 K km s , a very high emission level. When we use 
a 2 K km s _1 contour, clouds MC-6, MC-7 and MC-8 
in field 4C become one complex with a total molecular 
mass of 4.3xlO 5 Af . Similarly, in field 4B clouds 11 and 
12 (M tot =2.4xl0 5 M o ) and in field 2C clouds 16 and 17 
(M tot =2.1xlO 5 M0) become single complexes. Also, all 
the clouds (except cloud 3 which is at a very different ve- 
locity) in Field 1A merge into one complex. With better 
sensitivity, clouds 1, 2, 4 and 5 might be seen as belong- 
ing to one complex which would have a total mass of at 
least 4.9xl0 5 M Q . Although the difference in the mass 
of the merged clouds is not significantly different, the 
earlier conclusion that no cloud larger than 4xl0 5 M Q 
exists in M33 appears to be due to a different choice 
in cloud identification method and the low sensitivity of 
the observations. Our analysis is intended to show the 
need for using the same cloud identification and analysis 
technique for comparing cloud properties. 

In summary, though we only have a handful of clouds 
to compare with the Milky Way, we find that the M31 
and M33 GMCs are similar in most respects to those 
found in the Milky Way. The small number of clouds 
prevents us from studying the mass function of GMCs in 
M31. 



4.2. Automated Clump Finding Algorithms 

5 In fact, when we used Clumpfind on the Wilson & Scoville 
data, we found every one of their clouds, but Clumpfind actually 
merged 4 pairs of these clouds into a single cloud. 



Automated algo rithm s such as Clumpfind 
(Willi ams et al.l Il994h and Gaussclump 

( Stutzki fc GustenI I1990D work on the opposite phi- 
losophy from the single contouring algorithm described 
above: both Clumpfind and Gaussclump break apart 
a complex into as many parts as possible. Clumpfind 
contours the cube into discrete intervals specified by 
the user. Then it works its way through the contouring 
levels identifying isolated peaks as possible individual 
clumps. Clumpfind then assigns all pixels above the 
lowest contouring interval to one or the other clump. 
Finally it merges clumps which are closer than the 
spatial or velocity resolution. The advantage of this 
algorithm is that it has no preconceptions of the shape 
of the clouds and hence readily deals with the complex 
structure of GMCs. On the other hand, the contouring 
technique tends to clip away too many pixels which can 
lead to incorrect calculation of sizes and linewidths. 
This may be particularly problematic for data cubes 
with low dynamic range. Gaussclump avoids such 
problems of noisy or low dynamic range data by directly 
fitting clouds in the data cube. It tries to fit peaks of 
emission with an 11-parameter gaussian, also taking into 
account the resolution of the observations. The biggest 
disadvantage of this method is that it is restricted to 
defining clumps as gaussian, which is seldom the shape 
of the observed molecular emission. 

Neither method is particularly well-suited to analyse 
GMC properties. The reason has more to do with molec- 
ular clouds than the algorithms themselves. Molecular 
cloud structure is extremely complex, at least as observed 
in the Milky Way. These algorithms tend to divide a 
complex into as many subdivisions as allowed by the res- 
olution of the data. Hence these methods seldom find 
structures which are larger than than one or two resolu- 
tion elements. For example, if the resolution of a survey 
is 5 pc, these methods will not find a 100 pc complex 
because peaks of emission located more than 5 pc apart 
will be classified as separate clouds. The bias thus in- 
troduced is not fatal because cloud/clump properties in 
equivalent datasets (i.e. datasets with similar resolutions 
and noise characteristics) may still be compared, but gen- 
eralizations about a galaxy's cloud population can not 
be drawn. This is especially pertinent to extragalactic 
astronomy where higher and higher resolution observa- 
tions are continually being achieved with improvements 
in millimeter wave interferometry. 

We advocate that since molecular clouds generally have 
fairly sharp edges, the simplest cloud identification al- 
gorithms is the integrated intensity contouring method 
where one identifies the contour level associated with the 
edge of the cloud (e.g., S anders et al.l 119851 : iDame et al.l 
1986) in a 3 dimensional cube. The obvious shortcom- 
ing of this method is that it ignores the structure in the 
clouds. Another obvious disadvantage of the method is 
that it is difficult to automate and is quite time consum- 
ing because one has to make three different moment maps 
along the different axes and decide what the distinct 
complexes are. This method may suffer from blended 
emission bu t this confusi o n will be limited in extragalac- 
tic sources. iRosolowskyl (|2006h has recently introduced 
a new method for identifying and analyzing molecular 
clouds which may avoid many of the biases inherent 
in other methods discussed above; this method may be 



() 



Sheth et al. 



(J 

1 

i 


R I* 


o 

Si, 

1 

J** 




□ 


t s e , 30 12 

FWH« (Arcsocs) 




1 WhM (A'csacs) 


| 


























+ GEM 0B1 






ii y * * 








i 








05 2 







Fig. 4. — The three panels show the three derived properties 
compared with each other for each of the datasets: M31, M33 and 
the simulated Milky Way clouds, Cas A, Gem OBI and W3CO. 
No attempt has been made to deconvolve the lincwidths or radii 
because all datasets have the same linewidths and approximately 
the same angular resolution. The sizes are in arcseconds, linewidths 
in km s — 1 and amplitude in Jy beam -1 . There is no significant 
difference in the properties of the clouds identified by Gaussclumps. 
The largest M33 clouds are the least affected by the coarse (~ 
9" Resolution. The cluster of small M31 clouds with low linewidths 
are those in Complex C. These are very weak clouds and their 
properties are not well-defined. 

a useful tool for future surveys of molecular clouds in 
nearby spirals with CARMA and ALMA. 

4.3. M31, M33 and Milky Way Cloud Properties using 
Gaussclump 

We ran Gaussclump on all CO datasets in all three 
galaxies. We find 17 clouds in our BIMA M31 s urvey, 35 
clouds in the M33 survey of IWilson & Scovilld (|1990f ). 7 
clouds in the simulated Gem OBI complex, 5 clouds in 
the W3CO complex, and 13 clouds in the Cas A complex. 
These results should be compared to the six complexes 
found in M31 and the two complexes found in each of 
the simulated Milky Way clouds using the integrated in- 
tensity contouring method. The mean values for the am- 
plitudes, velocity widths and sizes are listed in Table El 
The velocity widths and sizes have not been deconvolved. 
Note that the M33 clouds were not simulated at the dis- 
tance of M31 and observed with BIMA. These data are 
similar to M31 but many fields have a slightly coarser 
resolution (8-9") than the BIMA data. 

The derived properties arc also shown graphically in 
Figure [4j The plots show no significant differences in the 
amplitudes, sizes or velocity widths for the three galaxies. 
We plot each of the derived properties against another 
derived property to test whether there is any dependence 
of one on the other, and find none. 

Our conclusion from this experiment is similar to that 
derived from the analysis of the integrated intensity con- 
touring method. The molecular clouds in these datasets, 
when compared in a consistent manner, show no signifi- 
cant differences. 

5. CONCLUSIONS 

In the northeastern spiral arm of M31, we have iden- 
tified six distinct, large complexes of molecular gas. All 



of these complexes lie along the spiral arm dust lanes, or 
in one case (Complex E), along a dust spur in between 
the dust lanes. These complexes have a mean diameter 
of 57±13 pc, a mean velocity width of 6.5±1.2 km s -1 , 
and a mean molecular mass of 3. Oil. 6 x 10 5 M©and 
are i ndistinguishable fro m those found in the Milky Way 
e.g- JPame et al.l (|1986ft . 

Meaningful comparison of GMCs in different galaxies 
requires consistent analysis of data taken with differ- 
ent instruments and different cloud identification tech- 
niques. This paper represents an attempt at eliminat- 
ing such differences and comparing GMCs in M31, M33 
and the Milky Way in as consistent a manner as possi- 
ble. We have simulated three Milky Way complexes at 
the distance of M31 and observed them with the BIMA 
array. The simulations show that interferometers are ex- 
cellent at recovering molecular complexes in the Local 
Group galaxies. We compared the simulated Milky Way 
complexes to the M31 data using an integrated inten- 
sity contour method and found that the complexes thus 
identified and analyzed f ell on the same lin ewidth size 
relationship as found by (|Pame et al.lll986l ). The M33 
clouds if analyzed in this manner also yield larger clouds 
than previously stated. Larger interferometric surveys of 
these galaxies are necessary to compare the GMC distri- 
bution and mass funct i on to that in the M ilky Way (e.g., 
lEngargiola eTa!][200l fcosolowskvlE2007th 

Finally, we compared two automated algorithms, 
Clumpfind and Gaussclump, and found that Gaussclump 
was better at recovering cloud properties than Clumpfind 
for low dynamic range data. Our experiments revealed 
that such algorithms are resolution dependent and the 
inherent clumpy nature of molecular emission prevents 
these methods from identifying clouds larger than one or 
two resolution elements. Hence we caution against us- 
ing these methods to characterize GMC populations in 
galaxies. Still, these may be used to quickly compare 
properties of clouds in data with similar noise and res- 
olution characteristics. Using Gaussclump, we find that 
the cloud properties in M31, M33 and simulated Milky 
Way data are indistinguishable. 



We are grateful for the insightful and helpful com- 
ments by the referee Elly Berkhuijsen which greatly im- 
proved the manuscript. We thank Leo Blitz, Greg En- 
gargiola, Tamara Heifer, Andrew Harris, Lee Mundy, 
Marc Pound, Erik Rosolowsky, Peter Teuben, Stephen 
White and Friedrich Wyrowski for helpful discussions. 
We acknowledge the use of NASA's Sky View facility 
(http://skyview.gsfc.nasa.gov) located at NASA God- 
dard Space Flight Center. We also acknowledge support 
for the BIMA array from the National Science Founda- 
tion under grant AST-9981289 to U. Maryland. 

Appendix 



Mode 



GMCs in M31,M33 & Milky Way 

Image After BIMA Model - BIMA 











m * 










1 1 1 i 1 1 1 1 i 1 'i 

1 1 ,11:1 ,l 



0.5 



1 



1 .5 



2 




Fig. 5. — The top row of this figure shows a channel from the simulation of the Gem OBI complex and the bottom row shows a channel 
from the simulation of the Carina complex. The first column is the Milky Way cloud at the distance of M31. The second column is the cloud 
simulated with interferometric observations with BIMA. The grey scale wedge below the rows is an indication of the flux in Jy beam - for 
the first three columns; it is a ratio for the final column. In column 2, one can immediately notice that the Carina cloud is significantly 
dimmer than its input model whereas the Gem OBI cloud appears unchanged. The third column is the model cloud minus the simulated 
image, and the fourth column is the simulated image divided by the model image. Column 3 is in Jy beam -1 . Notice that almost no flux 
is lost in the Gem OBI simulation whereas Carina arm cloud is weaker by as much as 30-50 %. 



APPENDIX 

SIMULATION OF BIMA OBSERVATION OF MILKY WAY GMCS AT THE DISTANCE OF M31 

The first step is to simulate the appearance of the Milky Way clouds if they were located at the distance of 
Andromeda. In other words, we need to scale the Milky Way maps to the distance of M31, which we did by scaling 
the map header values for the angular size of the map pixels by the ratio of the distance to each Milky Way complex 
to that of M31. The next step is to simulate observation of the distance-scaled Milky Way data with the BIMA 
interferometer. To get a precise match in the point-spread function, we modified the header values for the RA and 
Dec at the field center so as to match the position of the BIMA M31 observation; this enables the uv tracks for the 
simulated observations to match the actual observations. The original temperature scale which was in the units of T^ 
was converted to Jy beam -1 and the maps were gridded to a 2" pixel scale. A separate model was then created for 
each mosaicked pointing; each model had a primary beam taper applied to it. All seven of these models were then 
observed with the uv tracks used for M31. Since the original Milky Way data were approximately half as noisy as the 
BIMA data, we simulated additive noise in the uv plane and added that to our final maps. The final simulated Milky 
Way maps are therefore reflective of the effects of spatial filtering by the interferometer, and are similar in almost 
every way to the M31 and M33 data. 

EFFECT OF INTEROMETRIC FILTERING ON SIMULATIONS OF MILKY WAY OBSERVATIONS 

We initially began this comparative study project with a section of the Carina arm from the 1.2m survey because, 
at a heliocentric distance of 13.6 kpc, the Milky Way survey's linear resolution matched the BIMA observations of 
M31. However when we simulated the Carina clouds, the simulated maps showed significant differences as shown in 
Figure [5) The bottom row of this figure shows a single velocity channel from the Carina simulations. The first column 
is the model Carina complex at the distance of M31 before BIMA observations. The second column is the model after 
BIMA observations, i.e. the simulated map. Both of these columns are on the same greyscale and have the same 
contour levels. Column 3 is the input model minus the simulated map, and column 4 is the simulated map divided by 
the input model. One immediately sees the reduction in the flux in these maps. The clouds are dimmer by a factor 
of 30-50%. The flux that is resolved out is not a uniform background but rather structure on all scales. We further 
tested these effects by changing the orientation of the Carina model clouds by shifting them in different directions and 
rotating them at various angles. In every case, we saw the same result: the flux was reduced by as much as 50%. 

When we repeated the experiments for the Cas A, W3CO and Gem OBI clouds from the same Milky Way Survey, 
we found a completely different result. The interferometer recovers almost all of the flux for these simulations. An 
example of this is shown in the top row, which shows a single velocity channel from the Gem OBI simulations. The 
difference between the Carina clouds and these complexes is that these complexes are 4-6 times closer than Carina 
and were observed with a linear resolution of ~5 pc compared to a resolution of ^30 pc. The 1.2m maps of the Carina 
clouds therefore have smoother emission on a larger scale than the Cas A, W3CO and Gem OBI clouds; this smooth 
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TABLE 1 

Observing Parameters for M31 



a(J2000) a 


h 45 m 07 s .57 


<5(J2000) a 


41°35'46^'68 


Vlsr( km s _1 ) 


-135 


Total bandwidth (km s —1 ) 


142.85 


Velocity resolution (km s ) 


0.254 b 


Projected baselines (kA 


2.2-30.2 


Single-sideband Tgys (Kelvin) 


275 - 1000 


Calibration quasar 


0136+478 


Absolute flux calibrator 


Uranus 


Typical uncertainty in flux calibration 


15% 


Beam 


7'.' 11 x 5'.' 71 x 4.7° 


Noise in 2.6 km s — ^xhannel (Jy beam -1 ) 


0.15 



a Central pointing 

k Correlator data taken at 0.254 km s _1 resolution. Data were later 
smoothed to 2.6 km s -1 resolution for imaging. 

emission is then subsequently resolved out by t he interferometer. This filtering of flux f rom large, smooth structures by 
interferometers is expected and observed (e.g.. lBeck Sz Hoerneslll996tlHelfer et al.lr2002D . Presumably, the Carina arm 
at higher resolution would look like the other Milky Way complexes we studied (Gem OBI, W3CO and Cas A). Since 
the Carina arm model itself is incorrect, we cannot use the simulated Carina clouds in our comparative analysis. The 
recovery of the Gem OBI, W3CO and Cas A clouds is reassuring because it means that interferometers are excellent 
at recovering the molecular emission in GMCs; this is due to the inherently clumpy nature of GMCs. 
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TABLE 2 

Datasets & Instruments 



Galaxy 


Observatory 


Angular/Linear Resolution 


Assumed distance 


M31 


BIMA 


T.'ll X 5'.' 71 / 19-24 pc 


690 kpc 


M33 


OVRO 


5-7" / 19-26 pc 


790 kpc 


GEM OBI (Milky Way) 


CfA 1.2m 


8.8' / 5.1 pc 


2 kpc 


Cas A (Milky Way) 


CfA 1.2m 


8.8' 5.1 pc / 


2 kpc 


W3 (Milky Way) 


CfA 1.2m 


8.8' / 6.4 pc 


2.5 kpc 


Carina 3 - (Milky Way) 


CfA 1.2m 


8.8' / 34.8 pc 


13.6 kpc 



a Carina data not used in comparative study 



TABLE 3 
Properties of M31 GMCs 



Cloud 


RA 


DEC 


Vpcak 


V fwhm ' 


rp* C 

1 B 


Sco a 


Diameter 3 ' d 


M VIR ^ c 


M i a ' f 

ivi mol 




(J2000) 


(J2000) 


km s~ 1 


km s _1 


K 


Jy km s 


pc 


xl0 5 M© 


xl0 5 M Q 


A 


0:45:12.42 


41:37:07.23 


-126 


7.8 


3.32 


40.8 


44.5 


2.68 


3.13 


B 


0:45:12.73 


41:35:24.49 


-139 


7.8 


3.61 


67.6 


71.2 


4.29 


5.18 


C 


0:45:08.69 


41:27:49.58 


-134 




1.57 










D 


0:45:06.04 


41:35:23.25 


-134 


4.7 


2.02 


11.1 


41.1 


0.89 


0.85 


E 


0:45:04.97 


41:36:11.69 


-126 


6.2 


3.34 


55.4 


73.6 


2.80 


4.25 


F 


0:45:01.96 


41:36:38.99 


-131 


6.2 


4.24 


23.6 


57.4 


2.2 


1.80 


3 Typical errors 


arc: D,± 4pc, V f w hm>± 1-3 km s 1 , 


Sco, ±20%, Mvir,±30%, M mol 


,±30% 







Full width at half maximum measured from the spectrum 
c Brightest pixel in channel maps 

Average of deconvolved major and minor axis gaussian fit to integrated intensity map of complex 
6 Mvjh= 99 x V 2 f whm x Diameter 
f M mo i = 1.61 x 10 4 x d 2 x 

Sco This assumes a Galactic value for the conversion factor between the CO flux and the hydrogen 
column density of 3x 10 20 cm -2 (K km s -1 ) -1 and includes a factor of 1.36 to account for helium and other heavy elements. 



TABLE 4 

Properties of Milky Way GMCs Simulated at M31 



Cloud 


V fwhm 


rp* C 

L B 


Sco a 


Diameter 3 ' d 


My,/' c 


M mo; 3 ' f 




km s~ 1 


K 


Jy km s 1 


pc 


xl0 5 M Q 


xlO 5 M 


Geml 


7.2 


4.0 


70.87 


66.0 


3.4 


5.4 


Gem2 


11.1 


2.9 


105.4 


98.2 


12.0 


8.0 


Wl 


9.9 


2.4 


42.62 


40.1 


3.9 


3.3 


W2 


5.5 


1.3 


22.21 


61.5 


1.8 


1.7 


CasAl 


7.7 


4.5 


260.9 


123.0 


7.2 


20 


CasA2 


10.6 


2.9 


81.8 


62.5 


6.9 


6.3 



3 Typical errors are: D,± 4pc, V fwhm ,± 1.3 km s -1 , S C O,±20%, Mv/K,±30%, M mo! ,±30% 
k Full width at half maximum measured from the spectrum 
c Brightest pixel in channel maps 

d Average of deconvolved major and minor axis gaussian fit to integrated intensity map of 
complex 

° Mvir= 99 x V 2 fmhm x Diameter 
f M moi = 1.61 x 10 4 x d 2 x Sco 



TABLE 5 

Mean of Gaussclump Derived Properties 



Galaxy 


# of clouds 


Amplitude 


^ fwhm 


Size 






Jy beam -1 


km s _1 


Arcseconds 


M31 


17 


0.9±0.4 


3.6±1.0 


7.4±1.3 


M33 


35 


0.9±0.3 


4.9±1.7 


8.4±1.4 


Gem OBI 3 


7 


1.0±0.4 


4.3±0.7 


8.4±1.1 


W3CO a 


5 


0.6±0.2 


4.7±1.3 


7.6±1.3 


Cas A 3 


13 


1.0±0.4 


5.3±1.5 


8.8±1.8 



a These values refer to the simulated Milky Way clouds at the distance 
of M31. 



